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We study the electronic structure and correlations of vitamin B 12 (cyanocobalamine) by using the framework 
of the multi-orbital single-impurity Haldane-Anderson model of a transition-metal impurity in a semiconductor 
host. The parameters of the effective Haldane-Anderson model are obtained within the Hartree-Fock (HF) 
approximation. The quantum Monte Carlo (QMC) technique is then used to calculate the one-electron and 
magnetic correlation functions of this effective model. We observe that new states form inside the semiconductor 
gap found by HF due to the intra-orbital Coulomb interaction at the impurity 3d orbitals. In particular, the 
lowest unoccupied states correspond to an impurity bound state, which consists of states from mainly the CN 
axial ligand and the corring ring as well as the Co ^^-like orbitals. We also observe that the Co{3d) orbitals 
can develop antiferromagnetic correlations with the surrounding atoms depending on the filling of the impurity 
bound states. In addition, we make comparisons of the HF-i-QMC data with the density functional theory 
calculations. We also discuss the photoabsorption spectrum of cyanocobalamine. 


I. INTRODUCTION 

A. Metalloproteins and metalloenzymes 

Proteins are important building blocks of organisms. About 
one third of proteins have a metal atom such as iron, cobalt 
or zinc, and also about half of enzymes contain a metal 
atom^“^. The metal atom in these organometallic molecules is 
usually coordinated with nitrogen, oxygen or sulphur atoms. 
This kind of proteins and enzymes are called as metallopro¬ 
teins and metalloenzymes. An example for metalloproteins is 
hemoglobin which contains iron, and is involved in the trans¬ 
portation of oxygen and carbon dioxide in blood. An example 
for metalloenzymes is vitamin B12, which contains a cobalt 
atom and plays an important role in the production of red 
blood cells and in the functioning of the nervous system and 
the brain. 

Metalloenzymes and metalloproteins form an unusual class 
of organometallic molecules; they enable various chemical re¬ 
actions to take place in the biological environment of the or¬ 
ganisms instead of requiring, for example, higher values for 
pressure or temperature. The mechanisms for the function¬ 
ing of metalloenzymes and metalloproteins have been studied 
for more than half a century. However, much still remains to 
be understood about these organometallic molecules. In par¬ 
ticular, we think that it is important to know whether their 
unique chemical functions have a common underlying elec¬ 
tronic mechanism. 

As an example for metalloenzymes, in this paper we study 
vitamin B12, which is also called cobalamin (Cbl). We have 
chosen to study Cbl because it contains relatively few atoms 
and is suitable for numerical calculations in addition to its bi¬ 
ological importance. The molecular structure of Cbl was de¬ 
termined by X-ray measurements^, and consists of a ligand 
and a nucleotide attached to a corrin ring as shown in Fig. 
1 (a). In the corrin ring, four nitrogen atoms are attached to the 
cobalt atom which is located at the center. Here, we concen¬ 
trate on the cofactor cyanocobalamin (CNCbl) with the chem¬ 


ical formula C63H88C0N14O14P, in which a CN molecule is 
attached to Co as the axial ligand. The other two cofactors of 
Cbl are adenosylcobalamine (AdoCbl) and methylcobalamin 
(MeCbl), in which either an adenosyl or a methyl CH3 group 
is attached to Co instead of CN as the axial ligand. 


B. Comparison with the dilute magnetic semiconductors 

The CNCbl exhibits an energy gap of about 2.2 eV in its 
electronic spectrum^ and contains a transition metal atom. In 
this respect, it is similar to an entirely different class of com¬ 
pounds, the dilute magnetic semiconductors (DMS)^’^. The 
DMS materials such as (Ga,Mn)As are obtained by substi¬ 
tuting transition metal impurities into a semiconducting host 
material. They exhibit interesting physical properties such as 
long-range ferromagnetic correlations, which are found to ex¬ 
ist between the Mn magnetic moments at elevated tempera¬ 
tures. The transport measurements have shown that, in the 
dilute limit, an impurity bound state exists near the top of 
the valence band in the semiconductor gap^. The impurity 
bound state, which is a sharp resonant state in the single¬ 
particle spectrum consisting of spectral weight from the Mn 
impurity and the host, plays a crucial role in determining the 
electronic and magnetic properties of (Ga,Mn)As^^“^^ when 
studied within the framework of the Haldane-Anderson^^ 
model. For example, the long-range ferromagnetic correla¬ 
tions among the Mn impurities disappear rapidly as the impu¬ 
rity bound state becomes occupied with electrons. 


C. Motivation 

The comparison with the DMS materials leads naturally to 
the following questions: Do similar many-body effects play a 
role in determining the electronic and magnetic properties of 
CNCbl as well as the functioning of metalloenzymes and met¬ 
alloproteins in general? What is the special role of the tran- 
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FIG. 1. (Color online) (a) Schematic plot of the molecular structure 
of cyanocobalamine (C 63 H 88 C 0 N 14 O 14 P). (b) Schematic plot of the 
simpler truncated structure for cyanocobalamine, Im-[Co^^\corrin)]- 
CN^, for which we perform numerical calculations. 


sition metal atom in the organometallic molecules? In met- 
alloenzymes and metalloproteins, are there electronic states 
which are analogous to the impurity bound state found in the 
DMS materials? If so, what role to they play in the catalytic 
functioning of these molecules? The research presented in this 
paper is motivated mainly by these questions. 


D. Haldane-Anderson model 

As for the DMS materials, we study the electronic structure 
and correlations of CNCbl by using the Haldane-Anderson 
model of a transition metal impurity in a semiconducting 
host. For this purpose, we first map the electronic struc¬ 
ture of CNCbl onto the multi-orbital single-impurity Haldane- 
Anderson model by making use of the Hartree-Fock (HF) 
approximation. Then we study this effective Anderson 
Hamiltonian^^ by using QMC calculations with the Hirsch- 
Fye algorithm^^. Hence, we study CNCbl by combining the 
HF and the QMC techniques, which is called as HF-i-QMC. 
There have been many previous studies of the electronic state 
of Cbl and other metalloproteins by using the density func¬ 
tional theory Here, we study Cbl from the perspective 
of many-body physics within the framework of the Haldane- 
Anderson impurity model by using the HF-i-QMC approach. 

The multi-orbital single-impurity Anderson Hamiltonian^"^ 


is given by 

m,<7 y,G 

^yo^ma) ( 1 ) 

m,y,G 

V 

where (CmG) creates (annihilates) an electron in host state 
m with spin a, dl (dy) is the creation (annihilation) opera¬ 
tor for a localized electron at the Co(3dy) orbital, and riya = 
dlcydya. Here, £rn and e^y are the energies of the host and the 
3dy impurity states, respectively, Vmy is the hybridization ma¬ 
trix element between these states and Uy is the intra-orbital 
Coulomb repulsion. We have introduced a chemical poten¬ 
tial jl since the QMC calculations are performed in the grand 
canonical ensemble. 

This form of the Hamiltonian does not include the inter¬ 
orbital Coulomb interaction along with the Hund’s coupling. 
In cobaltates, the Hund’s coupling plays an important role^^, 
and we expect it to be important in Cbl also. In the present 
HF-fQMC calculations the inter-orbital Coulomb correlations 
are taken into account only at the HF level, and not in the 
QMC calculations. 

We note that the charge and spin states of Fe in hemoglobin 
were studied previously by using mean-field theory within the 
extended Haldane-Anderson model which contains the inter¬ 
orbital Coulomb interaction and the Hund’s coupling^^. In 
addition, the high-spin to low-spin transition due to oxygen 
binding in myglobin was studied by using mean-field theory 
for the extended Haldane-Anderson model^^. 


E. Truncated molecular structure for cyanocobalamine 

We perform these calculations by using the simpli¬ 
fied molecular structure Im-[Co^^^(corrin)]-CN+ shown in 
Fig. 1(b), as has been done in previous DFT calculations^"^. 
Here, it is clearly seen that Co is coordinated by 5 N and 1 
C atoms, which forms an approximately octahedral config¬ 
uration. In order to obtain this simplified structure, the nu¬ 
cleotide group attached to the corrin ring is truncated and the 
side chains of the corrin ring are replaced by hydrogen atoms. 
The resulting Im-[Co^^^(corrin)]-CN+ has 56 atoms, instead of 
181 in CNCbl, and contains 238 electrons instead of 718. Our 
use of this truncated structure for modelling CNCbl decreases 
the number of the host eigenstates and the computational cost 
of the QMC calculations, because in the HF-fQMC approach 
it is necessary to calculate all of the host Green’s functions in 
order to determine the chemical potential. 


F. Outline of the paper 

The outline of this paper is as follows. In Section II, in 
order to obtain the one-electron parameters of the Anderson 
Hamiltonian, we perform Hartree-Fock calculations with the 
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Gaussian program^ In particular, we first write the HF so¬ 
lution for the Fock matrix in the orthogonal basis set of the 
natural atomic orbitals (NAO)^^’^^. From this, we obtain the 
one-electron parameters e^y, and Vmv^ The HF calculation 
finds an energy gap of 8.5 eV for the host electronic states. 
In addition, the hybridization matrix elements are largest be¬ 
tween the Co ^g-like states and the atoms around Co, in par¬ 
ticular, the CN axial ligand. For the intra-orbital Coulomb re¬ 
pulsion f/y, we use the value obtained from the two-electron 
integrals. However, we also perform calculations for different 
values t/y in order to see the dependence of the results. In 
the HF-fQMC approach, the intra-orbital Coulomb repulsion 
is taken into account by both the HF and the QMC parts of 
the calculations. In order to eliminate this double-counting, 
the bare Co(3Jy) energy levels are shifted by a constant 
amount. 

In Section III, we present QMC data on the single-particle 
and magnetic correlation functions for the effective Ander¬ 
son Hamiltonian. The Hartree-Fock calculations find that 
the host eigenstates have a semiconducting gap in the en¬ 
ergy interval — 10.2eV < £ < —1.7 eV. The HF-fQMC re¬ 
sults presented in Section III show that new single-particle 
states are induced in the intervals —lOeV^e^ —8.5eV and 
—5.5 eV^ e ^ — 2 eV. We find that the states in the inter¬ 
val — lOeV^e^ — 8.5 eV arise mainly from the doubly- 
occupied Co r 2 g-like states, while containing smaller amount 
of states from the corrin ring. On the other hand, based 
on the QMC data on the single-particle and magnetic corre¬ 
lation functions, we conclude that the states in the interval 
—5.5eV^e^ —2eV represent impurity bound states arising 
from the strong hybridization of the Co ^^-like states with the 
host eigenstates from the valence band. These impurity bound 
states have sufficient spectral weight to accommodate 2 elec¬ 
trons including spin. Interestingly, only 0.4 of these electrons 
originate from the Co ^^-like orbitals, the remaining 1.6 elec¬ 
trons originate from the host states which consist of orbitals 
mainly from the CN axial ligand and to a lesser extent from the 
corrin ring. In Section III, we also present QMC data on the 
magnetic correlation functions. We study the magnetic mo¬ 
ment formation as a function of the electron filling. In addi¬ 
tion, we show that magnetic moments develop in the valence- 
band host eigenstates which are strongly hybridized with the 
Co ^g-like orbitals. The size of these magnetic moments de¬ 
pend sensitively on the filling of the impurity bound states 
located at — 5.5 eV and « —4 eV. In addition, the magnetic 
moments of these host states are coupled antiferromagneti- 
cally to the Co ^^-like magnetic moments, and these correla¬ 
tions vanish rapidly with the electron filling of the impurity 
bound states. It is because of this filling dependence that we 
attribute the induced states at « — 5.5 eV and ^ —4 eV to be 
arising from impurity bound states. It would be useful to ex¬ 
perimentally look for signatures of the antiferromagnetic cor¬ 
relations which are predicted by the HF-fQMC calculations. 
It would also be useful to look for experimental evidence that 
the states in the interval — 5.5 eV ^ ^ —2 eV indeed arise 

from impurity bound states in CNCbl. 

In Section IV, we present a discussion in which we compare 
the HF-fQMC data with the DFT calculations on the same 


molecule. With respect to the HF results, the DFT calcula¬ 
tions find that new single-particle spectral weight is induced 
in the intervals —10.5 eV < £ < —8.0 eV and —5.5 eV < £ < 
— 1.7 eV. This is similar to the HF-fQMC data. However, 
the overall distributions of the spectral weight obtained by the 
DFT and the HF-i-QMC calculations are very different. In the 
DFT calculations, the states in the interval —10.5 eV < £ < 
—8.0eV do not arise from the upper-Hubbard states of the 
doubly-occupied Co t 2 g orbitals. In addition, the states in the 
interval —5.5eV < £ < —1.7eV do not correspond to an im¬ 
purity bound state. Furthermore, the magnetic correlations are 
absent in the DFT calculations. In Section IV, we also com¬ 
pare the HF-fQMC data with the experimental data on the pho¬ 
toabsorption spectrum of CNCbl, which is characterized by 
several distinctive peaks. The HF-fQMC data suggest that the 
lowest-energy excitations are dominated by the electron trans¬ 
fer from the Co r 2 g-like states to the impurity bound states. 
In Section IV, we also discuss that the combined DFT-fQMC 
technique may also be applied to the same problem. In addi¬ 
tion, we note that the inter-orbital Coulomb interaction with 
the Hund’s coupling needs to be included in the QMC calcu¬ 
lations before carrying out quantitative comparisons with the 
experimental data on CNCbl. 

Section V gives the summary and conclusions of the paper. 


II. CONSTRUCTION OF THE EFFECTIVE ANDERSON 
HAMILTONIAN 

Various methods for the ah initio calculation of the pa¬ 
rameters of the Anderson Hamiltonian had been developed 
for rare-earth, actinide and 3J compounds^"^. Here, we use 
the Hartree-Fock approximation to determine the parameters 
of the effective Haldane-Anderson model for the truncated 
molecule Im-[Co^^^(corrin)]-CN+. 

In Section II. A, we discuss how the self-consistent Hartree- 
Fock solution for the Fock matrix is obtained. In Section II.B, 
we introduce the natural atomic orbitals (NAO’s) which form 
an orthogonal basis^^. In Section II.C, we use the Fock matrix 
written in the NAO basis to determine the bare energy lev¬ 
els and the hybridization parameters of the effective Haldane- 
Anderson model. In Section II.D, we present results on the 
density of states obtained from the Hartree-Fock approxima¬ 
tion. By mapping to the effective Haldane-Anderson model, 
the energies of the host states £rn and the impurity levels e^v 
are obtained. Here, we show the density of states of the host 
part of the model along with the positions of the bare Co(3 J) 
levels. The hybridization matrix elements obtained from the 
mapping to the Anderson Hamiltonian are presented in Sec¬ 
tion ILF. In Section ILF, we discuss the values of the onsite 
Coulomb repulsion used in the QMC simulations. In order to 
eliminate the effects of the double counting of the Coulomb 
interaction in HF-fQMC, the Co(3 Jy) energy levels are shifted 
by /iDC. which is discussed in Section II.G. The results pre¬ 
sented in Section II are based on Ref. [35]. 

At this point, it is necessary to discuss how the atomic co¬ 
ordinates of the molecule Im-[Co^^^(corrin)]-CN+, which are 
needed as input parameters for the Hartree-Fock calculations. 
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are determined. Here, we use the approach taken by Kornobis 
et which may be summarized as follows. X-ray mea¬ 
surements had been used to obtain with high accuracy the ge¬ 
ometrical parameters for the molecule MeCbl^^. In order to 
have an initial estimate of the atomic coordinates of CNCbl, 
the axial ligand CH 3 in MeCbl was replaced by CN. The re¬ 
sulting molecular structure for CNCbl was then truncated ob¬ 
taining the simplified structure Im-[Co^^^(corrin)]-CN+. The 
final atomic coordinates were determined after performing an 
optimization on the coordinates of the simplified structure 
by using the DFT technique with the 6-31G(d) Gaussian ba¬ 
sis set^^ and the BP 86 exchange functionals^^’^^. This way 
we obtain a reasonable set of values for the geometrical pa¬ 
rameters, which are used as input in the Hartree-Fock calcu¬ 
lations. The atomic coordinates used in this paper for Im- 
[Co^^^(corrin)]-CN+ are the same as in Ref. [24]. 


A. Hartree-Fock solution and the Fock matrix 


According to the electronic structure theory, the many-body 
wave function 'F of the molecule is given by the Slater deter¬ 
minant of the one-electron molecular orbitals y/n. Within the 
Hartree-Fock mean-field approximation^^, the one-electron 
wave function is determined from the self-consistent 

solution of 



|ym(rOP 

|r-r'| 


Vn{r) 


( 2 ) 


y f j3 I 

Ir-r'| 




F are given by 


pij = 


(3) 


where 


N,i 

CyijCfi j 


( 6 ) 


forms the density matrix. Here, Nq\ denotes the total num¬ 
ber of the electrons and also twice the number of the occu¬ 
pied molecular orbitals. For Im-[Co^^^(corrin)]-CN+, we have 
Nq\ = 238. Using as input the one- and two-electron integrals, 
the solution of the Roothan equation for and the density 
matrix are obtained self-consistently through iteration. In or¬ 
der to obtain the HF solution for Im-[Co^^^(corrin)]-CN+, we 
have used the Gaussian program^ ^ with the 6-3IG basis set 
and N = 347 basis functions. We perform the HF calculations 
with the Gaussian 6-3IG basis set instead of 6-31G(d), which 
was used in the initial optimization, in order to increase the 
speed of the QMC calculations reported in Section III. Af¬ 
ter the HF solution has been obtained, we have formed the 
Fock matrix. From the Fock matrix we have obtained the 
one-electron parameters of the effective Anderson Hamilto¬ 
nian through a procedure outlined in the following two sub¬ 
sections. 


B. Basis of the natural atomic orbitals 


with En the molecular-orbital energy. Here, the left-hand 
side includes the kinetic-energy and the electron-ion potential 
Fion(r) terms in addition to the Hartree and the Fock contribu¬ 
tions. 

The molecular orbitals can be expressed in terms of 
the N atomic orbitals |0;), 

N 

\\lfn)='£Cni\(l)i), (3) 

i 

where Cni are the elements of the coefficient matrix C. By 
substituting Eq. (3) into Eq. (2), the Roothan equation is 
obtained^^, 

CF = ECS. (4) 

Here, the overlap matrix S has the matrix elements Sij = 
{(t>i\(t>j)- Since the atomic orbitals |0/) do not form an orthog¬ 
onal basis, the overlap matrix is not diagonal. The matrix E is 
diagonal with the diagonal elements giving the molecular or¬ 
bital energies, Enm = SnmEm- The elements of the Eock matrix 


The overlap matrix S defined above is not diagonal, since 
the atomic orbitals do not form an orthogonal basis set. A 
convenient orthogonal basis for expressing is provided 
by the basis of the natural atomic orbitals The NAO’s 

form a maximally localized basis set; they are similar to the 
atomic orbitals but they are also orthogonal. The NAO’s are 
obtained by minimizing the function 

fd^r\^i{r)-^i{r)\^, (7) 

i 

where Wi denotes the electron occupation number of the 
atomic orbital \(j)i) obtained by the Hartree-Eock calculation. 
Here, we will first express the Eock matrix in the NAO basis, 
and from this deduce the parameters of the Anderson Hamil¬ 
tonian except for the Coulomb repulsion. 

Using the basis of the NAO’s, the molecular orbitals \\i/n) 
can be expressed as 

5 N-5 

\¥n)='£Pnv\^dv)+'£Pni\ti), ( 8 ) 

v=l i=l 

where, |fev)’s represent the five NAO’s for the Co(3dy) or¬ 
bitals, and |0i)’s represent the remaining NAO’s. We have 
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obtained the elements of the Fock matrix within the NAO ba¬ 
sis by using the Gaussian program^ ^ along with the NBO 6.0 
software package^^. 


C. One-electron parameters of the Anderson Hamiltonian 
from the Fock matrix 

Within the HF-fQMC approach, we determine the one- 
electron parameters of the Anderson Hamiltonian from the 
matrix elements of the Fock matrix expressed in the orthog¬ 
onal NAO basis. This procedure is sketched in Fig. 2(a) and 
(b). Here, the A x A Fock matrix F consists of the 3 J part 
F 3 J, the host part Fq, and the hybridization parts M and its 
transpose as illustrated in Fig. 2(a). The host part of the 
Fock matrix Fq is not diagonal in the NAO basis. By diagonal¬ 
izing Fq, we obtain the host energy levels of the effective 
Haldane-Anderson model, and the corresponding host eigen¬ 
states \um)- In turn, the hybridization matrix elements Vmv be¬ 
tween the host eigenstates lum) and the Co(3dy) NAO’s \^dy) 
are obtained. We denote by F' the Fock matrix expressed in 
the basis formed from the Co(3dy) NAO’s and the orthogonal 
host states, which is illustrated in Fig. 2(b). We determine the 
one-electron parameters of the Anderson Hamiltonian from 
F'. The impurity energies e^v of the effective model are given 
by the diagonal matrix elements of F 3 J. We neglect the off- 
diagonal elements of F 3 J, which are denoted by A,v' Fig- 
2 (b), since they are negligible in comparison to the diagonal 
terms. 

Using the procedure outlined above, we obtain the model 
parameters and Vmv from the transformed Fock ma¬ 

trix F'. Considering the effective Hamiltonian Eq. (1), we 
note that now (dya) creates (annihilates) a fermion in the 
Co(3dy) NAO state |0jv), while (Cma) creates (annihi¬ 
lates) a fermion in the orthogonal host eigenstate \um)- Hence, 
the Vmv’s represent the hybridization matrix elements between 
the host eigenstates \um) and the Co(3dy) NAO’s. 

There have been previous studies where the electronic 
structure calculations are combined with the many-body 
techniques'^^. The use of NAO’s here for the case of a 
molecule is similar to the use of maximally-localized Wan- 
nier orbitals'^^ to obtain the model parameters for various 
compounds'^^’'^'^. 
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FIG. 2. (a) Sketch of the Fock matrix F in the natural atomic or¬ 

bital basis. Here, the 3d and the host parts of the Fock matrix are 
separated as F^d and Fq, respectively. The hybridization parts are 
M and M^. The host part Fq is not diagonal in the NAO basis, (b) 
Sketch of the Fock matrix F' which is obtained after diagonalizing 
the host part. We obtain the one-electron energies and the hybridiza¬ 
tion matrix elements of the effective Haldane-Anderson model from 
F'. Eigenvalues of the host part of the Fock matrix F' yields the 
host energy levels of the Anderson model. The off-diagonal matrix 
elements yield the hybridization parameters Vmy between the host 
states and the Co(3dy) NAO’s. Since tyyf are negligible, e^/y’s be¬ 
come the energy levels of the impurity states. 


D. Host density of states and the Co{Sdy) energy levels 

Figure 3(a) shows the total density of states D{£) defined 
by 

D{e) = f^5{e-En) (9) 

n=l 

where En are the eigenvalues of the Fock matrix. Here, 
the highest occupied molecular orbital (HOMO) is located 
at - 10.2 eV, and the lowest unoccupied molecular orbital 
(LUMO) is at -1.7 eV. Hence, the occupied and the unoccu¬ 
pied states are separated by an energy gap of 8.5 eV. In Fig. 


3(b), the density of the host eigenstates defined by 

N-5 

( 10 ) 

m=l 

where £m are the energy levels of the host eigenstates, is 
shown. We observe that in Fig. 3(b) the HOMO and the 
LUMO levels are slightly shifted in comparison to Fig. 3(a). 
The energy levels of the Co(3Jy) NAO states, ejy, are also 
shown as vertical lines in Fig. 3(b). We observe that the xz, 
— y^, and yz orbitals are degenerate and are much lower in 
energy with respect to the 3z^ — and xy orbitals, which are 
located above the energy gap and near the LUMO level. 
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E. Host-Co(3dv) hybridization 


In this section, we present Hartree-Fock data on the hy¬ 
bridization matrix elements of the Haldane-Anderson model, 
which were obtained as discussed in Subsection ILC. The 
square of the hybridization matrix elements IVmvP between 
the m’th host eigenstate and the Co(3dy) NAO are shown 
as a function of the host energy £fn in Fig. 4. Here, we observe 
that the magnitude of IVmv P the cobalt 3z^ — and xy or¬ 
bitals can be an order of magnitude larger than those for xz, 
x^ — and yz. This is not surprising since the 3z^ — and xy 
orbitals lie along the Co-N and Co-C bond directions accord¬ 
ing to the coordinate system used here. Figure 4(a) shows 
that, for the V = 3z^ — orbital, IVmvP has the largest val¬ 
ues at m = 114 and 116. For v = xy, the largest values occur 
at m = 111 and 112. In Section III, we will see that the host 
eigenstates m = 114, 111, 112 and 116 are the states which 
are most affected by the presence of the cobalt atom. 

We note that IVmv P’s have large values for a few host states, 
while being negligibly small for the remainder. This feature 
may be used in the future to increase the speed of the numeri¬ 
cal calculations without losing accuracy. 


F. Coulomb repulsion at tbe Co(3dv) orbitals 


FIG. 3. (Color online) (a) Total density of states of Im- 
[Co^^\corrin)]-CN+, D{e), obtained within the Hartree-Fock ap¬ 
proximation by using the Gaussian program. The vertical solid and 
dashed lines denote the HOMO and LUMO levels, respectively, (b) 
Density of states of the host eigenstates of the effective Haldane- 
Anderson model, D/j(e). Here, the HOMO and LUMO levels are 
slightly different than in (a). In addition, the energy levels of the 
Co(3dy) NAO’s are indicated by the vertical lines. 


In CNCbl, the 5 nitrogen atoms and the C atom surround¬ 
ing Co form roughly a local octahedral environment for Co. 
Hence, the Co(3dy) levels separate into two groups, one con¬ 
taining the 3z^ — and xy orbitals and the other containing the 
xz, x^ — y^, and yz orbitals. This is similar to the separation of 
the 3d orbitals into the groups of Cg and t 2 g orbitals in a crys¬ 
tal field with cubic symmetry. However, here we have taken 
a coordinate system such that the x and y axis are located at 
45 degrees to the Co-N bond direction in the corrin layer in¬ 
stead of being parallel to the Co-N bond direction. Hence, the 
x^ — y^ and xy orbitals have been exchanged. In the follow¬ 
ing, because of our choice of the coordinate system which is 
the same as that of the Gaussian program, we will denote the 
3z^ — and xy orbitals as the ^^-like states, and the xz, x^ — y^, 
and yz orbitals as the t 2 g-likQ states. 


The intra-orbital Coulomb repulsion at the Co(3dy) atomic 
orbitals is given by the two-electron integral 

t/v = y d^r j |fev(rOP- (H) 

Using the 6-3IG Gaussian basis set^p we obtain Uy = 
36.8 eV for the Co(3dy) orbitals. However, we note that the 
value of Uy depends on which Gaussian basis set is used. 
While for the 3-21 basis set Uy = 42.8 eV, for the 6-31G and 
the 6-31G(d) sets, we obtain Uy = 36.8 eV. 

Since the effective Anderson Hamiltonian introduced above 
uses the basis of the Co(3dy) NAO’s, the Coulomb two- 
electron integrals also need to be evaluated in this basis in¬ 
stead of the basis of the Co(3dy) atomic orbitals. When we 
evaluate the inter-orbital Coulomb repulsion given by Eq. (11) 
using the Co(3dy) NAO’s with the Gaussian 6-3IG basis set, 
we obtain « 28 eV. 

In the remaining of this paper, we will use a constant U for 
Uy, and present QMC results particularly for U varying be¬ 
tween 28 eV and 36 eV. We are neglecting the inter-orbital 
Coulomb interactions in the QMC calculations, which are ac¬ 
tually close in value to the intra-orbital terms. Hence, in the 
QMC part of the calculations, we are taking into account the 
effects of the Coulomb repulsion at a simple level, treating U 
more as a variable than a true ab initio parameter. Inclusion 
of the inter-orbital Coulomb interaction will clearly provide 
more accurate results as discussed in Subsection IV.D. 
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FIG. 4. (Color online) Hartree-Fock results on the square of the hy¬ 
bridization matrix elements | P between the host eigenstates |w^) 
and the Co(3Jv) NAO states \^dv) plotted as a function of the energy 
of the m’th host eigenstate In (a) results are shown for the eg- 
like V = 3z^ — and xy orbitals, and in (b) for the t 2 ^-like v = xz, 
x^ and yz. Here, the vertical solid and dashed lines denote the 
HOMO and LUMO levels, respectively. We observe that the host 
states m= 114, 111, 112 and 116 have the strongest hybridization 
matrix elements. The matrix elements with the t 2 ^-like states are 
much smaller than those with the ^^-like states. 
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FIG. 5. (Color online) Density of states of the host states of the effec¬ 
tive Haldane-Anderson model D/j(e) from the Hartree-Fock calcula¬ 
tion. Here, the Co(3 Jy) levels have been shifted by jioc in order to 
compensate for the double-counting of the Coulomb repulsion by the 
HF-i-QMC approach. The shift has been evaluated for U = 36 
eV. 


It is possible to argue that, because we are using only 
the intra-orbital Coulomb repulsion in QMC, the double¬ 
counting shift should be orbital dependent and given by = 
U (^j^)/2, where ) is the orbital dependent Co(3Jy) oc¬ 
cupation number including spin from the HF calculation. 
However, if we include the inter-orbital Coulomb interaction 
in QMC, then the shift would be orbital independent. Since 
we try to model the actual system as much as possible, we use 
a constant /Idc- 

In Fig. 5, the shifted Co(3Jy) NAO energy levels are 
shown for f/ = 36 eV and /ioc = 0.68f/. Also, shown here 
is the host density of states Dh{£). For these parameters, the 
xz, x'^ — y^ and yz states are located below the bottom of the 
valence band, while the 3z^ — and xy states are located in 
the middle. 


G. Double counting in HF-hQMC 

At this point, it is necessary to note that in the HF-fQMC ap¬ 
proach the onsite Coulomb repulsion U is taken into account 
twice, once in the Hartree-Fock calculation, and a second time 
in the QMC calculations. In order to prevent the double¬ 
counting of the effects of U in the many-body approaches 
which combine the electronic-structure calculations with the 
QMC simulations, a constant /Idc is usually subtracted from 
the Co(3Jy) levels"^^"^^, £dv = ^dv — Mdc- Here, /Idc is 

given by 

MDC = t/^, (12) 

where {n^^) is the total number of the electrons (includ¬ 
ing spin) in the Co(3d) levels obtained within Hartree-Fock. 
Since the Hartree-Fock calculations yield (n^^) = 6.8, we use 
Mdc = 0.681/ in this paper. 


III. QUANTUM MONTE CARLO RESULTS 

In this section, we present QMC data on the effec¬ 
tive Haldane-Anderson model for Im-[Co^^^(corrin)]-CN+, of 
which parameters were obtained in the previous section. For 
this model, we have performed QMC calculations by using 
the Hirsch-Fye QMC algorithm^^. In particular, we have cal¬ 
culated the Co(3 Jy) and the host single-particle Green’s func¬ 
tions and the magnetic correlations functions. 

In Subsection III.A, we introduce the correlation functions 
which we evaluate with QMC. In Subsection III.B, we present 
data on the electron occupation number and the magnetic mo¬ 
ment formation at the Co(3Jy) NAO’s. Here, we observe that 
in-gap states develop arising from the Co(3Jy) states in the 
semiconductor energy gap. The induced in-gap states reduce 
the value of the semiconductor gap found by the Hartree-Fock 
calculation. We find that the Co(3Jy) states induced above 
the top of the valence band correspond to the upper Hub- 
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bard states of the Co r 2 g-like orbitals. On the other hand, the 
Co(3 Jv) states induced in the middle of the semiconductor gap 
correspond to the impurity bound states originating from the 
Co eg-likc states. We reach this assignment based on the QMC 
results on the chemical potential dependence of the Co(3dy) 
electron occupation number and magnetic moment formation, 
host electron occupation numbers and magnetic moment for¬ 
mation as well as the antiferromagnetic correlations between 
the Co eg states and the host eigenstates. 

The QMC results on the host electrons are shown in Sub¬ 
section III.C. Here, we observe that the host electrons form 
new states inside the semiconductor gap at the same ener¬ 
gies where the Co(3dy) NAO’s develop in-gap states. These 
host states are the ones which have the strongest hybridization 
matrix elements with the Co ^^-like states. These host states 
also develop magnetic moments which depend strongly on the 
electron filling of the impurity bound states. 

In Section III.D, we present data on the antiferromagnetic 
correlations which develop around the Co atom. We find that 
the antiferromagnetic correlations develop mainly between 
the Co eg-like states and the electron spins at the CN axial 
ligand, and they depend strongly on the filling of the impurity 
bound states. The results presented in Section III are based on 
Ref. [51]. 

The QMC calculations are performed within the grand 
canonical ensemble, and the total electron number {nr) is ob¬ 
tained by evaluating with QMC, 

(ht) = ii(^ va^va) + IK (13) 

v=l <7 m=l <y 

In presenting the QMC data in the next section, we will obtain 
results as a function of fi and evaluate the corresponding {nr) 
also with QMC. For the truncated structure Im-[Co^^^(corrin)]- 
CN+, the total electron number is Nq\ = 238. By comparing 
{ht) with A/ei, we obtain the corresponding value for fi. 

A. QMC measurements 

The Matsubara-time dependent single-particle Green’s 
function for the Co(3dy) NAO states is defined by 

Gy^{t) = -{T,dya{t)dl^{0)), (14) 

where is the usual Matsubara T-ordering operator and 
dyai'^) = e^^dy(je~^^. Similarly, the host Green’s function 
is defined by 

G„<t(t) = -(7Vc„<,(t)4<,(0)). (15) 

In addition, we calculate the square of the magnetic moment 
at each Co(3dy) orbital, ((M^)^), where 

My = 4t^vt - ■ (16) 

In order to probe the magnetic correlations around the Co 
atom, we also calculate the equal-time Co(3 Jv)-host magnetic 
correlation function where 

— ^lrv\^rn-\ ~ ( 1 ^) 


We present data on as well. 

In the following subsections, we show QMC data on the 
Co(3dy) electron number {uy) and the square of the Co(3dy) 
moments ((M^)^) as a function of the chemical potential /i. 
We also present data on the total host electron number {rih) 
as well as the total electron number {nr) as a function of 
fl. In addition, we study the local moment formation at the 
host eigenstates and their magnetic coupling to the Co(3dy) 
NAO’s. We find that antiferromagnetic correlations develop 
between the Co Cg states and the electronic spins at the CN 
axial ligand depending on the value of fl. 

In obtaining the QMC data shown below, Matsubara-time 
discrete steps between At = 0.04 eV“^ and 0.07 eV“^ were 
used. The results are presented for temperature T = 2000^ in 
the grand canonical ensemble. Hence, a temperature broaden¬ 
ing of order 0.17 eV is expected for the energy resolution of 
the QMC data. 


B. QMC results on the Co(3dv) electrons 

We begin presenting QMC data with Figure 6(a) which 
shows the electron number {uy) of the Co(3dy) NAO states 
as a function of the chemical potential jl for U = 36 eV. Sim¬ 
ilarly, Figure 6(b) shows the square of the magnetic moment 
((M^)^) at the Co(3dy) orbitals as a function of fi. We observe 
that the results for the Co(3d) xz, — y^, and yz states are sim¬ 
ilar to each other. The results for the 3z^ — and xy states are 
also similar. From the fl dependence shown in these figures it 
is possible to extract information on the single-particle spec¬ 
tral weight distribution and the local moment formation. 

The bare energy levels e^v of the Co(3d) xz, x^ — y^, and yz 
NAO states are almost degenerate and are located at e ^ — 45 
eV below the lower edge of the valence band, as seen in 
Fig. 5. In Fig. 6(a) we see that for these orbitals {uy) ex¬ 
hibits a jump of unit magnitude at /i « —45 eV. These orbitals 
remain singly occupied as /i is increased upto the interval 
— lOeV^/i^ —8 eV, at which point these orbitals become 
nearly doubly occupied. In Fig. 6(b), we see that the mag¬ 
netic moments at these orbitals have their maximum value in 
the interval —45 eV^ /i ^ — 10 eV. Above « —10 eV, these 
moments decrease rapidly due to double occupancy. The in¬ 
gap states in the energy interval —lOeV^e^ —8eV have 
sufficient spectral weight to accommodate a total of nearly 
three electrons, and their location corresponds to approxi¬ 
mately ejv +1/. Hence, we deduce that these in-gap states 
correspond to the upper-Hubbard states of the Co t 2 g-likQ or¬ 
bitals. We note that they are located inside the gap, and hence 
new impurity states are formed inside the gap found by HF. 

Figure 5 shows that the bare energy levels £dv of the 3z^ — 
and xy orbitals are near —25 eV. In Fig. 6(a) we observe that 
for these two orbitals, (riy) becomes finite at /i ^ — 38 eV, 
increases rapidly around fl ^ —30 eV, and remains less than 
one until fi reaches —5.5 eV. This type of ji dependence 
suggests that the f/ = 0 single-particle spectral weight, which 
is a 5-function near —25 eV, has been shifted down in en¬ 
ergy and has been broadened. These orbitals remain less than 
singly occupied until fl reaches ^ — 5.5 eV, where (rixy) ex- 
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FIG. 6. (Color online) (a) QMC results on the electron occupation 
number {riy) of the Co(3 Jy) natural atomic orbitals plotted as a func¬ 
tion of the chemical potential /i. (b) Square of the magnetic moment 
{{My)^) at the Co(3dy) natural atomic orbitals versus the chemical 
potential /i. Here, the vertical solid and dashed lines denote the 
HOMO and the LUMO levels, respectively. These results are for 
U = 36 eV. 


hibits a sudden increase of about 0.2. We think that this sud¬ 
den increase corresponds to an impurity bound state as found 
in the mean-field solution of the Haldane-Anderson model^^. 
For the 3z^ — orbital, a similar increase occurs at a slightly 
higher energy of ^ —4 eV. We note that these jumps take 
place in the energy gap found by HF. Hence, new impurity 
states are formed in the energy gap. As seen in Fig. 6(b), the 
magnetic moments at the xy and 3z^ — orbitals also exhibit 
sudden increases at these energies. The occupation numbers 
of these states remain constant at unit magnitude as jl is in¬ 
creased upto « 10 eV. This value of fl corresponds approxi¬ 
mately to £dv + U for the xy and 3z^ — orbitals. Above this 
energy, the xy and 3z^ — orbitals rapidly become doubly 
occupied, and the magnetic moments decrease. 

In summary, from Fig. 6(a) we deduce that new in-gap 
states originating from the Co(3J) xz, and yz or¬ 

bitals form in the energy interval —lOeV^e^ —8 eV, while 
new states from the xy and 3z^ — orbitals form at energies 
~ — 5.5 eV and ~ —4 eV, respectively. These newly formed 
states reduce the value of the semiconductor gap found by the 
Hartree-Fock calculation, which is 8.5 eV, down to ^3 eV. 

Figure 7(a) shows the QMC data on the total electron occu- 
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M (eV) 

FIG. 7. (Color online) (a) Total electron occupation number (n^) 
of the Co(3d) natural atomic orbitals as a function of the chemical 
potential fi. (b) Total number of the host electrons (rih) versus jU. 
(c) Total number of electrons, {nr) = {n^) + (w/z), versus fd. Here, 
results are shown for various values of the intra-orbital Coulomb re¬ 
pulsion U. In addition, the vertical solid and dashed lines denote the 
Hartree-Fock results for the HOMO and LUMO levels, respectively. 

pation of the Co(3dy) states 

5 

M = L Zi^ladvc) (18) 

V=1 o 

versus the chemical potential /i. Here, results are shown for 
various values of U in order to see the dependence on U. We 
note that in these calculations, as U is varied, also varies 
according to Eq. (12). 

For U = 36 eV, we observe that {ud) increases as /i is var¬ 
ied from -10 eV to -8 eV, because the xz, x^ — y^ and yz orbitals 
become nearly doubly occupied in this interval. In the energy 
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gap, for —8 eV^ /i ^ — 5.5 eV, we get {rid) ~ 7.55. In com¬ 
parison, the X-ray scattering experiments find that the value 
of (rid) in CNCbl is 7.7 ± 0.1^^. Above ^ -5.5 eV, we ob¬ 
serve two step-like increases in (ud) coming from the xy and 
3z^ — states. For U = 36 eV, the energy gap between the 
highest occupied and the lowest unoccupied Co(3 J) levels is 
« 3 eV. In comparison, the experimental value for the semi¬ 
conductor gap of CNCbl is about 2.2 eV^. As U is decreased 
from 36 eV to 28 eV, we observe that the in-gap states move 
closer to the semiconductor gap edges, and the magnitude of 
the energy gap increases. 

Also shown in Fig. 7(a) are results for U = 0, in which 
case we obtain = 6.8 in the interval —12eV</i<0 eV. 
As seen in Fig. 3(b), for U = 0, the bare energy level for the 
xz, x^ — and yz orbitals is ^ — 20 eV, while for the 3z^ — 
and xy orbitals one has « — 1 eV and « 0 eV, respectively. 
If the hybridization is turned on while keeping U = 0, then 
spectral weight of the 3z^ — and xy orbitals are broadened 
and shifted to higher energies. For this reason, (rid) remains 
constant in the interval —12 eV < /i < 0 eV for f/ = 0. 

Finally, we note that the hybridization matrix elements 
y^v’s are dependent on the atomic coordinates used in the HF 
calculation. Because of uncertainties in the atomic coordi¬ 
nates, it is necessary to study how the general features of the 
single-particle spectrum depend on V^y. With this in mind, 
we have repeated the above QMC calculations by reducing 
the hybridization matrix elements by 10%. In this case, we 
found that, for U = 36 eV, the lowest energy impurity bound 
state occurs at —6.5 eV instead of —5.5 eV, and the energy 
gap is reduced to ^ 2 eV. Hence, we see that the impurity 
bound state is robust with respect to small variations in V^^y. 
This is important because the values of the depend on 
the Co-C bond length, and we are using simple estimates for 
the geometrical parameters as described in the introduction of 
Section II. 


C. QMC results on the host electrons 

Figure 7(b) shows the total number of the host electrons 

N-5 

{nh) = £ £(4cTCm(T), (19) 

m=l o 

where N = 347, plotted as a function of the chemical potential 
fl. Comparing Figs. 7(a) and (b), we observe that (rih) and 
(rid) exhibit jumps at the same values of fi. Hence, the host 
states also develop spectral weight inside the gap at the same 
energies as the Co(3Jy) states. We note that the error bars 
in the data on (uh) are larger than in (^j), since the QMC 
calculation of (rih) is computationally more costly. 

For the interval —lOeV^e^ —8 eV, the jump in (ud) is 
« 2.5, while the jump in (^/j) is « 0.5. On the other hand, 
we note that, in the interval —5.5eV^e^ —2 eV, the jump 
in (rih) is about 1.6, while that in (rid) is « 0.4. Hence, 
actually there are more host states induced in the gap for « 
—5.5 eV ^ e ^ — 2 eV than from the Co(3Jy) orbitals. 



FIG. 8. (Color Online) HF+QMC results on the energy gap Eg be¬ 
tween the HOMO and LUMO states versus the intra-orbital Coulomb 
repulsion U (left axis). Total electron number at the Co(3d) orbitals, 
{rid) when the chemical potential is in the energy gap, versus U (right 
axis). Here, ^Udc which shifts the bare Co(3Jy) energy levels is var¬ 
ied along with U according to Eq. (12). 


The total number of the electrons {nr) = {rid) + {^h) ver¬ 
sus fl is shown in Fig. 7(c). For Im-[Co^^^(corrin)]-CN+, 
the total electron number is 238. From this figure we iden¬ 
tify /i ^ — 8.5 eV as the new value of the HOMO level, and 
ji ^ —5.5 eV as the new LUMO level. Here, we also observe 
that the in-gap states located in the interval —5.5eV^/i^ — 

3.5 eV have sufficient spectral weight to accommodate a to¬ 
tal of two electrons. In comparison, the in-gap states in 
— lOeV^/i^ — 8 eV accommodate 3 electrons. 

We also note that £dy for 3z^ — and xy states are near —25 
eV. Hence, for U = 36 eV, the doubly-occupied states occur 
at /i « ejy + U = 11 eV. Hence, the in-gap states which form 
from the 3z^ — and xy orbitals do not originate from the 
upper Hubbard states. Instead, they are similar in origin to the 
impurity bound states found in the mean-field solution of the 
Haldane-Anderson model of a transition-metal impurity in a 
semiconductor host^^. 

Since for Im-[Co^^^(corrin)]-CN+ the total electron number 
is 238, in Fig. 7(c) we see that the effective energy gap de¬ 
creases as U is increased. For U = 36 eV, the effective en¬ 
ergy gap becomes « 3 eV. Hence, the energy gap can be sig¬ 
nificantly reduced with respect to the Hartree-Fock value of 

8.5 eV. 

The values of the effective energy gap are plotted in Fig. 
8 as a function of U. We note that in this calculation, as U is 
varied, the shift /ioc is also varied according to Eq. (12). Also 
plotted in Fig. 8 is {rid) versus U. For {rid), experimen¬ 
tal estimate^^ is 7.7 ±0.1. Hence, for U ^36 eV, HF-fQMC 
yields values for the energy gap and {rid) which are compara¬ 
ble to the experimental values. For this reason, in the remain¬ 
der of this paper we will present QMC data for U = 36 eV. 

We next explore the nature of the in-gap states originating 
from the host eigenstates. For this purpose, in Figure 9(a) we 
show the occupation of the host eigenstates {rim) as a func¬ 
tion of the host energy 8^ for fi = — 7.0 eV and U = 36 eV. 



11 


2 

-1-1- 

1 — ^ —n—^ — I 


112« t 

.116 1 - 


• 

1 n=-7 eV ■ 


m=114 

1 



1 

5 1 


1 

1 

1 

n 

^^ 

1 1 

w 

- • 




1 ' n ' 

0.6 


i (b) : 


m=114* 

1 

1 

0.4 

r 

1 

1 

1 

0.2 

112* ' 

1 

'll6 1 

1 


• 

1 

n 


1^ - 


-20 -10 0 10 


(eV) 


FIG. 9. (Color Online) (a) Occupation number of the host eigenstates 
{rim) versus the host energy £m for chemical potential fi = —1.0 eV. 
Here, we see that the host eigenstates m= 111, 112, 114 and 116 are 
not doubly occupied even though the corresponding £m are located 
deep below jU. (b) Difference in the occupation number of the m’th 
host eigenstate as /i is changed from — 7.0 eV to — 2.2 eV, (Arim) = 
(^m)|g=-2.2eV “ (^m) |g=-7.oeV» plotted as a function of £m- These 
results are for U = 36 eV. In addition, here, the vertical solid and 
dashed lines denote the HOMO and the LUMO levels, respectively. 


Here, fl is located inside the energy gap found by QMC. In 
this figure, we notice that there are valence-band host eigen¬ 
states which are not doubly occupied even though they are 
located deep below the chemical potential. In particular, the 
host eigenstates with labels m = 114, 111, 112 and 116 ex¬ 
hibit the most reduced values of (um)- We note that these host 
eigenstates have the strongest hybridization matrix elements 
with the 3z^ — and xy orbitals, as it was seen in Fig. 4(a). 

Next, we increase ji from —7.0 eV to —2.2 eV, and plot 
in Fig. 9(b) the increase in (rim), (Arim) = (^m)|g=-2.2eV - 
(^m)|g=-7.0eV, as a function of the host energy £m in Fig. 
9(b). We already know from Fig. 7(b) that, for U = 36 eV, 
the in-gap states become fully occupied when /i = — 2.2 eV. 
Here, we see that the largest increase in {rim) takes place for 
m = 114, 111, 116 and 112. The rest of the host eigenstates 
exhibit negligible change as fl is increased from — 7.0 eV to 
—2.2 eV. Hence, the impurity bound states derive mainly from 
the m=114. 111, 116 and 112 host eigenstates. In turn, within 
the rigid band picture, it is possible to state that the host eigen¬ 
states m= 114, 111, 116 and 112 have significant amount of 
single-particle spectral weight between the energies — 7.0 eV 




FIG. 10. (Color online) (a) Top view of Im-[Co^^\corrin)]-CN+. We 
note that the x and y axes are located at 45 degrees to the Co-N bonds 
in the corrin plane, and the z axis is perpendicular to the corrin plane, 
(b) Side view of Im-[Co^^\corrin)]-CN^, where the imidazole and 
CN ligands attached to Co below and above the corrin ring, respec¬ 
tively, are seen. Here, the carbon and nitrogen atoms have been la¬ 
belled for convenience, while the hydrogen atoms are not shown. 


and —2.2 eV. 

In order to gain insight into the real-space structure of 
the in-gap host eigenstates located in the interval between 
—7.0 eV and —2.2 eV, we study the NAO composition of 
these host eigenstates. However, first, in Fig. 10 we dis¬ 
cuss the geometrical structure of Im-[Co^^^(corrin)]-CN+. Fig. 
10(a) shows a top view of the corrin ring around the Co atom 
in Im-[Co^^^(corrin)]-CN+. Here, we have only shown the C 
and N atoms around Co, and have not included the H atoms. 
According to our coordinate system, which is the same as 
that of the Gaussian program, the corrin ring lies in the xy 
plane, and the x and y axes make 45 degrees with the Co- 
N bond directions. Figure 10(b) shows the side view of Im- 
[Co^^^(corrin)]-CN+. Here, we see the imidazole ring attached 
to Co below the corrin plane and the CN axial ligand attached 
to Co above the corrin plane. Here, the C and N atoms in the 
axial ligand are labelled as 30C and 3IN as in the notation of 
the Gaussian program. 

Figure 11(a) shows the NAO composition of the m = 111’th 
host eigenstate. In particular, here we plot the NAO weight of 
the m = 111’th host eigenstate defined as 

\Umif = ( 20 ) 

versus the label i of the host NAO state. Figure 1 l(b)-(d) show 
similar results for the m = 112, 114 and 116’th host states. 
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FIG. 11. (Color online) NAO weight of the m’th host eigenstate 
\umi\^ = I (0/1 Wm) P versus the label i of the host NAO for (a) m = 111, 
(b) 112, (c) 114, and (d) 116. The maximum value for index i is 347, 
which is the total number of the host NAO’s. 


In these figures, we see that the m = 111 and 112’th states 
have large weight from the 2/7;r NAO’s of the CN axial ligand. 
Similarly, the m = 114 host state mainly consists of the 2pG 
NAO’s of the CN axial ligand. On the other hand, the m = 116 
host state consists of the 2pa NAO’s on the lO’th, 15’th and 
20’th C sites as well as N sites located around Co in the corrin 
plane as sketched in Fig. 10(a). Thus in Figs. 10 and 11, we 
see the real-space composition around the Co site of the in-gap 
host states induced in the interval —7.0 eV ^ /l ^ — 2.2 eV by 
the presence of Co. 

In Fig. 12, we show results on the occupation number and 



FIG. 12. (Color online) (a) Occupation of the m’th host eigenstate 
{rim) versus p. (b) Square of the magnetic moment of the m’th host 
eigenstate ((M^)^) versus p. Here, the vertical solid and dashed 
lines denote the HOMO and the LUMO levels, respectively. In addi¬ 
tion, these results are for U = 36 eV. 


the magnetic moment for the same host eigenstates m= 111, 
112, 114 and 116. Figure 12(a) shows the occupation num¬ 
ber (rim) versus p for these host eigenstates. We observe 
that a sharp increase in (rim) occurs at p ^ £m in the valance 
band. We see that the single-particle spectral weight of the 
m = lll’th host state is broadened in the valence band. In 
addition, a step-like increase exists inside the energy gap at 
« —5.5 eV, above which this host state reaches double occu¬ 
pancy. Similar p dependencies are observed for the others. 
The host state m = 112 also exhibits a jump at —5.5 eV, 
while for m = 114 and 116 the jump occurs at — 4 eV. The 
jumps in (rim) are largest for m = 114 and 111. 

Figure 12(b) shows the square of the magnetic moment 
for the same host eigenstates. We observe that sig¬ 
nificant size moments develop for these host states until the 
in-gap states become occupied. As the in-gap host states are 
filled and they approach double-occupancy, the magnetic mo¬ 
ments decrease rapidly. 

Finally, we note that the states induced in the interval 
— lOeV^e^ —8eV originate mainly from the Co t 2 g-likQ 
states. In Figs. 7(a)-(c), we observe that in this interval there 
is sufficient spectral weight to accommodate 3 electrons, of 
which about 2.5 electrons have the Co t 2 g-likQ character. The 
remaining 0.5 electrons originate from the m = 115 host state, 
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which mainly consists of the 2pK orbitals of the C and N 
atoms in the corrin ring. We do not show the composition 
of the m = 115’th host state in order not to increase the length 
of the manuscript. 


D. QMC results on the Co(3dv)-host magnetic correlations 

In Fig. 13, we discuss the magnetic correlations between 
the Co(3Jv) electrons and the host eigenstates m= 111, 112, 
114 and 116. In particular, these figures show the magnetic 
correlation function between the Co(3Jv) state and 

the m’th host eigenstate plotted as a function of /i. 

In Figs. 13(a) and (b) we observe that the Co{?>dxy) NAO 
develops antiferromagnetic correlations with the m = 111 and 
112’th host states. These antiferromagnetic correlations di¬ 
minish rapidly as the in-gap state located at ^ — 5.5 eV be¬ 
comes occupied. The m = 111 and 112’th states do not exhibit 
magnetic correlations with the other Co(3 Jy) orbitals. In Figs. 
13(c) and (d) we observe that the Co(3d^^2_p) NAO develops 
antiferromagnetic correlations with the m = 114 and 116’th 
host states. In Figs. 1 l(a)-(c), we have seen that the m = 111, 
112 and 114’th states contain large weight of the NAO’s local¬ 
ized at the CN axial ligand. Thus we deduce that the electrons 
at the CN axial ligand develop antiferromagnetic correlations 
with the Co Cg orbitals. These antiferromagnetic correlations 
diminish rapidly when the impurity bound states become filled 
with electrons. Hence, as Im-[Co^^^(corrin)]-CN+ goes from 
the Aei = 238 state to Aei = 239 and 240 states, the antifer¬ 
romagnetic correlations between the CN ligand and the Co Cg 
electrons disappear rapidly. We note that the m = 116 host 
state contains spectral weight from the 2pG NAO’s at the C 
and N sites in the corrin ring as seen in Fig. 11(d). Figure 
13(d) shows that these NAO’s also develop antiferromagnetic 
correlations with the Co(3J3^2_^2) orbital. 

It is important to note that we identify the in-gap induced 
states at ^ — 5.5 eV and ^ —4 eVas impurity bound states be¬ 
cause of the p dependencies discussed in Figs. 6, 7, 12 and 13. 
Similar behaviour is obtained in studying the impurity bound 
states found in the DMS materials 

We have also studied the magnetic correlations between the 
magnetic moments forming at the Co 3d^^2_p and 3dxy or¬ 
bitals. We find that these orbitals exhibit only weak magnetic 
correlations. This is different than the case for the DMS ma¬ 
terials, where an indirect ferromagnetic coupling is generated 
between two impurity moments due to their antiferromagnetic 
coupling to the same continuum of host electronsHere, 
the Co 3d^p_p and 3dxy orbitals do not have strong hybridiza¬ 
tion with the same host eigenstates as seen in Subsection II.E. 
In modelling the DMS materials, a continuum of host states 
with constant hybridization is used, while for CNCbl the host 
eigenstates have discrete energy levels with varying hybridiza¬ 
tion matrix elements. Hence, an indirect ferromagnetic cou¬ 
pling is not generated. However, the inclusion of the Hund’s 
coupling can generate a ferromagnetic coupling between the 
magnetic moments at the Co 3dy^2_p and 3dxy NAO’s. 

It would be interesting to probe experimentally the antifer¬ 
romagnetic correlations found in the HF-fQMC calculations. 



M (eV) 


FIG. 13. (Color online) Magnetic correlation function (MyM^) be¬ 
tween the m’th host eigenstate and the various the Co{3dy) NAO’s. 
Here, results are shown for host eigenstates (a) m = 111, (b) 112, 
(c) 114, and (d) 116. In addition, the vertical solid and dashed lines 
denote the HOMO and the LUMO levels, respectively. These results 
are for t/ = 36 eV. 


It is known that Cbl exhibits weak diamagnetism^^’^^, which 
has been intepreted as Co in Cbl having a low spin state. This 
is also partly responsible for the notion that the interaction ef¬ 
fects are not important in Cbl. In the HF-i-QMC results, we see 
that magnetic moments can develop in the Co Cg orbitals and 
in the host eigenstates coupled to these orbitals. It is possi¬ 
ble to calculate the magnetic susceptibility within HF-fQMC 
in order to extract the effective moment size for Co and the 
whole molecule. But such a calculation would be more mean¬ 
ingful after including the Hund’s coupling. 
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IV. DISCUSSION 

In this section, we discuss the meaning and implications 
of the HF+QMC data shown in Section III. For this pur¬ 
pose, in Subsection IV.A, we compare the HF-fQMC data 
with the DFT calculations on the same truncated molecule 
Im-[Co^^^(corrin)]-CN+. We find that there are important dif¬ 
ferences between the results of these two approaches. In 
particular, while the HF-i-QMC calculations for the Haldane- 
Anderson model find that the LUMO state corresponds to an 
impurity bound state, this is not the case for the DFT results. 
In Subsection IV.B, we compare the HF-i-QMC data with the 
photoabsorption experiments on CNCbl. Here, we discuss the 
nature of the lowest excited states in the photoabsorption spec¬ 
trum. The HF-fQMC data suggest that the lowest-energy exci¬ 
tations are dominated by electron transfer from the Co r 2 g-hke 
orbitals to the impurity bound states which contain spectral 
weight mainly from the CN axial ligand and the corrin ring. 
In Subsection IV.C, we note that the DFT-i-QMC method may 
also be applied to the same problem. In Subsection IV.D, 
we emphasize the importance of including the inter-orbital 
Coulomb repulsion and the Hund’s coupling. 


A. Comparison of the HF+QMC data with the DFT results 
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FIG. 14. (Color online) DFT results on the total density of states 
D(e) of Im-[Co^^\corrin)]-CN+ obtained by using the Gaussian pro¬ 
gram with the BP86 exchange potentials and the 6-3IG basis func¬ 
tions. The solid and dashed vertical red lines denote the HOMO 
and LUMO levels obtained by the DFT calculation. Similarly, the 
solid and dashed vertical black lines denote the HOMO and LUMO 
levels obtained by the Hartree-Fock calculation described in Sec¬ 
tion II. Here, we see that in the DFT calculations additional spec¬ 
tral weight is induced in the intervals —10.5 eV < e < —7.8 eV and 
-5.7 eV < e < —1.8 eV when compared with the HF results. 


Various DFT based calculations have been performed for 
Cbii6-26 subsection, we compare the HF+QMC data 

obtained for Im-[Co^^^ (corrin)]-CN+ by using the Haldane- 
Anderson model with the results obtained from the density 
functional theory (DFT) calculations on the same truncated 
molecule. We note that these DFT results are similar to the 
results of the time-dependent DFT calculations on the same 
system^"^. 

Figure 14 shows the total density of states D(e) versus £ 
for Im-[Co^^^(corrin)]-CN+ obtained by using the Gaussian 
program^^ with the 6-3IG basis functions and the BP86 ex¬ 
change potentials^ . Here, the solid and dashed vertical red 
lines denote the HOMO and LUMO levels obtained by DFT. 
According to the DFT results the energy gap is « 2.2 eV, 
which is in agreement with the experimental results. The 
vertical black lines, on the other hand, denote the HOMO 
and LUMO levels obtained by the HF calculations shown 
in Fig. 3(a). Comparing the HF results on D(e) shown 
in Fig. 3(a) with the DFT results of Fig. 14, we observe 
that the HF HOMO and LUMO levels located at —10.5 eV 
and — 1.7eV move to —7.8eV and —5.6eV, respectively. 
Hence, in comparison to the HF results, the DFT calcula¬ 
tion yields additional spectral weight in the energy intervals 
-10.5 eV<e< -7.8eVand-5.6eV<e< -1.7eV. 

Next, we compare the DFT results on D(e), Fig. 14, 
with the single-particle spectral weight deduced from the 
HF+QMC data on the total electron number {nr) versus ji 
shown in Fig. 7(c) for U = 36 eV. We observe that the loca¬ 
tions of the HOMO and LUMO levels obtained by DFT and 
HF+QMC are comparable. In addition, the DFT result for the 
total number of the electrons in the Co(3dy) NAO states is 7.4, 
which is similar to the HF+QMC result seen in Fig. 8. So, the 


DFT and the HF+QMC results for the energy gap and (rid) 
are comparable to each other as well as to the the experimen¬ 
tal results. However, we will see in the following that there 
are important differences between the DFT and HF+QMC re¬ 
sults when the overall distribution of the single-particle spec¬ 
tral weight near the semiconductor gap edges is considered. 

In the HF+QMC data presented in Section III for U = 36^y , 
we have seen that new states are induced in the energy in¬ 
terval — lOeV^e^ —8eV which originate mainly from the 
doubly-occupied Co 3dxz, 3d^2_y2 and 3dyz orbitals. This is 
because + U for the Co t 2 g-likQ orbitals corresponds to this 
energy interval. In the same interval, there also exist smaller 
amount of host states due to hybridization with the Co t 2 g- 
like states. On the other hand, the states induced in the in¬ 
terval —5.5 eV^ e ^ — 2 eV in the HF+QMC data originate 
from the Co 3dy^2_^2 and 3dxy states. These states are simi¬ 
lar in nature to the impurity bound states found in the mean- 
field and QMC studies of the Haldane-Anderson model for 
the DMS materials^^"^^. We have also seen that host states, 
which are mainly originating from the CN axial ligand, are 
induced in the energy interval — 5.5 eV ^ e ^ — 2 eV. In fact, 
here, the amount of the single-particle spectral weight origi¬ 
nating from the host states is larger than that from the ^^-like 
states. From Figs. 7(a) and (b), we observe that in the inter¬ 
val —5.5 eV ^ e ^ — 2 eV the host states accommodate about 
1.6 electrons while the ^^-like states accommodate about 0.4 
electrons. 

Now, we discuss the nature and origin of the states 
forming in the energy intervals —lOeV^e^ —8eV and 
—5 eV ^ e ^ — 2 eV in the DFT calculations. For this pur- 







15 


Molecular orbital 

n 

E„ (eV) 

3^2 _ ^2 

xy 

Co(3dv) 

xz 

2 2 

X — y 

yz 

30C 

(summed over 2p’s) 

31N 

(summed over 2p’s) 

HOMO-10 

109 

-10.22 

0.0 

0.0 

0.0 

0.0 

0.02 

0.01 

0.01 

HOMO-9 

110 

- 9.70 

0.0 

0.0 

0.10 

0.0 

0.01 

0.14 

0.14 

HOMO-8 

111 

- 9.69 

0.0 

0.0 

0.01 

0.0 

0.16 

0.25 

0.24 

HOMO-7 

112 

- 9.27 

0.04 

0.0 

0.0 

0.0 

0.0 

0.09 

0.49 

HOMO-6 

113 

- 9.19 

0.01 

0.0 

0.0 

0.0 

0.0 

0.14 

0.28 

HOMO-5 

114 

- 9.15 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.01 

HOMO-4 

115 

- 8.74 

0.0 

0.0 

0.0 

0.08 

0.15 

0.09 

0.21 

HOMO-3 

116 

- 8.43 

0.0 

0.02 

0.0 

0.83 

0.02 

0.01 

0.03 

HOMO-2 

117 

- 8.07 

0.0 

0.01 

0.66 

0.0 

0.0 

0.02 

0.12 

HOMO-1 

118 

- 7.80 

0.02 

0.0 

0.0 

0.0 

0.38 

0.01 

0.03 

HOMO 

119 

- 7.77 

0.01 

0.0 

0.01 

0.01 

0.10 

0.01 

0.02 

LUMO 

120 

- 5.64 

0.0 

0.0 

0.0 

0.0 

0.05 

0.0 

0.0 

LUMO+1 

121 

- 5.08 

0.03 

0.49 

0.02 

0.01 

0.0 

0.0 

0.0 

LUMO+2 

122 

- 4.87 

0.38 

0.07 

0.01 

0.0 

0.01 

0.06 

0.02 

LUMO+3 

123 

- 4.74 

0.09 

0.01 

0.04 

0.0 

0.0 

0.01 

0.01 

LUMO+4 

124 

- 4.07 

0.0 

0.01 

0.0 

0.0 

0.0 

0.0 

0.0 

LUMO+5 

125 

- 3.83 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

LUMO+6 

126 

- 2.93 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

LUMO+7 

127 

-2.20 

0.0 

0.0 

0.0 

0.0 

0.0 

0.01 

0.01 

LUMO+8 

128 

- 2.03 

0.0 

0.0 

0.0 

0.0 

0.0 

0.01 

0.01 

LUMO+9 

129 

- 1.93 

0.0 

0.0 

0.0 

0.0 

0.01 

0.01 

0.01 

LUMO+10 

130 

- 1.78 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


TABLE L DFT results on the NAO spectral weights |C„/p of the molecular orbitals near the semiconductor gap edges for Im-[Co^^\corrin)]- 
CN^. Here, n is the index for the molecular orbital, and i is for the NAO. The first two columns denote the molecular orbital and the 
corresponding index. The following five columns denote how much Co(3d) NAO spectral weight exists in the n'th molecular orbital. Similarly, 
the last two columns denote the amount of the total 2p spectral weight at the C and N atoms in the CN axial ligand. For convenience, when 
\Cni\'^ > 0 . 05 , the numbers are printed in bold face. 


pose, we present DFT results in Table I on the NAO spectral 
weights of the molecular orbitals located in these intervals. 
Using the DFT data, we expand the one-electron molecular 
orbitals \\f/n) in terms of the NAO’s of the whole molecule 
including both the host and the Co(3 J) states, 

N 

\Wn)=Y,C„i\^i)- ( 21 ) 

i 

Here, the sum over the NAO index i goes upto N=347. Ta¬ 
ble I shows DFT results on \Cni\'^ for the molecular orbitals 
which are located right below and above the DFT semicon¬ 
ductor gap. In particular, here results are shown for the molec¬ 
ular orbitals with index 109 < < 130, which includes the 

states located in the intervals —10.22 eV < e < —7.77 eV and 
-5.64 eV < £ < —1.78eV. We are interested in how much 
Co(3d) NAO weight these molecular orbitals contain. We are 
also interested in how much NAO weight of the CN axial lig¬ 
and is included in these energy intervals. We will compare 
them with the HF+QMC data. It is for these purposes that the 
data in Table I are presented. 

The first and the second columns of Table I denote the lo¬ 
cation of a given molecular orbital with respect to the semi¬ 
conductor gap and the molecular orbital index n. The third 


column shows the corresponding molecular energy En. The 
following five columns show \Cni\^ for the Co(3J) NAO states. 
The last two columns show the NAO weights summed over the 
three 2p orbitals, Y. 2 p IQ/P. foi* the C and N atoms of the CN 
axial ligand. We note that the C and N atoms of the CN axial 
ligand are labelled as 30C and 3 IN following the notation of 
the Gaussian program^ ^ 

Here, we observe that the molecular orbitals with the in¬ 
dex ^ = 115 through 118 contain significant amount of spec¬ 
tral weight from the Co 3dxz, ^d^ 2 _y 2 and 3dyz^ NAO’s. In 
comparison, the HOMO state {n= 119) and the LUMO state 
{n = 120) contain smaller amount of Co 3J^^ NAO weight. On 
the other hand, the molecular orbitals n = \2\ and 122 con¬ 
tain significant amount of Co 3dxy and 3d'^^2_^2 NAO weights, 
respectively. 

Next, we compare these DFT results with the HF+QMC 
data from Section III. In the DFT results, we see that in 
the interval —lOeV^e^ —8eV there is spectral weight 
originating from the Co 3d^2_y2, 3dxz and 3dyz NAO states. 
This agrees with the HF+QMC data. However, the amounts 
of the single-particle spectral weight contained in this inter¬ 
val are different. According to the DFT results, in this in¬ 
terval the Co 3dxz, 3d^2_y2 and 3dyz NAO’s accommodate 
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1.56, 1.84 and 1.88 electrons, respectively, including spin. 
On the other hand, as seen in Fig. 6(a), the HF+QMC ap¬ 
proach finds that each of these orbitals contains in this inter¬ 
val about 0.9 electrons including spin. In this case, the interval 
—10 eV^e ^ — 8 eV is approximately equal to the value of 
^dv + U. Hence, the states in this interval correspond mainly 
to the upper Hubbard states of the r 2 g-hke NAO’s. This com¬ 
parison shows that the DFT calculations do not yield the split¬ 
ting of the r 2 g-like states into the lower and upper Hubbard 
states found by HF-fQMC. 

Next, we discuss the DFT results for the interval 
—5eV^e^ —2 eV, which contains spectral weight from the 
Co 3djcy and 3dy^2_^2 NAO states. The HF-fQMC data also 
finds spectral weight from the Co 3dxy and 3dy^2_y.2 NAO’s 
in this interval. However, the amounts of the spectral weight 
contained in this interval are different. According to the DFT 
results, the Co 3dxy and 3dy^2_j.2 NAO’s contain 1.0 and 1.14 
electrons in this interval. The HF-fQMC yields 0.2 electrons 
for each of these orbitals in the same interval. 

In the last two columns of Table I on the DFT results, we 
observe that in the interval —5.5eV^e^ —2 eV, there is lit¬ 
tle amount of NAO spectral weight originating from the CN 
axial ligand. In the HF-fQMC data, however, there is sig¬ 
nificant amount of spectral weight from the CN axial lig¬ 
and in this energy interval. This is because the CN ax¬ 
ial ligand hybridizes strongly with the Co 3dxy and 3d^2_y2 
NAO’s. Hence, in the DFT calculations the states in the in¬ 
terval —5.5 eV^ e ^ — 2 eV do not arise from an impurity 
bound state as found in the HF-fQMC calculations. 

When we compare the DFT and the HF-i-QMC results, we 
find good agreement for the magnitude of the HOMO-LUMO 
gap and the total occupation number of the Co(3J) or¬ 
bitals. However, there exist important differences in the over¬ 
all distribution of the single-particle spectral weight near the 
semiconductor gap edges. In particular, we see that the DFT 
calculation yields much less spectral weight in the interval 
—5.6eV^e^—1.8eV originating from the CN axial ligand. 
In addition, according to DFT , the Co(3 Jy) spectral weight in 
the interval —lOeV <e<—8 eV does not arise from the split 
upper-Hubbard states of the Co r 2 g-like orbitals. 

Furthermore, in Table I we see that, in the interval 
— 10.22 eV < e < —1.78 eV, there are 22 molecular orbitals, 
which can accommodate a total of 44 electrons. On the other 
hand, in the HF-i-QMC data shown in Fig. 7(c), we observe 
that there is spectral weight to contain only 5 electrons in 
the interval —lOeV^e^ —2eV. Hence, the single-particle 
spectrum in this interval is much denser according to the DFT 
results. 

In addition to these, we have also performed DFT 
calculations by using the B3LYP exchange-correlation 
functional^^’^^ instead of the BP86. We do not present these 
results here because of the length of the manuscript. In this 
case, the energy gap is 3.2 eV. However, the distribution 
of the spectral weight near the HOMO and LUMO levels is 
similar to that shown in Table I for BP86. 

Furthermore, we have also performed calculations by us¬ 
ing the local spin density approximation^"^ (LSDA) for Im- 
[Co^^^(corrin)]-CN+ with the Gaussian program^ In that 


case, the energy gap is ^2.1 eV, and the overall distribution 
of the single-particle spectral weight is similar to the DFT re¬ 
sults shown in Table I. As in the DFT calculations, the LSDA 
approach does not capture the Hubbard splitting of the r 2 g-hke 
orbitals, and the states at the LUMO level do not arise from 
an impurity bound state. 

In addition to metalloproteins and metalloenzymes, the 
DFT technique is widely applied for studying transition- 
metal-containing organometallic molecules, which are used 
in fields ranging from solar-cell applications to hydrogen stor¬ 
age. The differences discussed here between the DFT and the 
HF-fQMC approaches might simply be due to the fact that 
the present HF-fQMC approach is too simple to describe the 
electronic state of CNCbl. If that is not the case, then these 
comparisons suggest that caution is necessary in interpret¬ 
ing the results of the DFT technique when used for studying 
organometallic molecules containing transition metal atoms. 

In addition to these, the DFT and HF-fQMC results differ 
on the nature of the electronic correlations. The HF-fQMC 
technique finds that antiferromagnetic correlations exist be¬ 
tween the Co eg magnetic moments and the electronic spins 
localized at the CN axial ligand depending on the filling of the 
impurity bound states. However, the magnetic correlations are 
not contained in the DFT results, as it is well known. 


B. Comparison of the HF+QMC data with the 
photoabsorption spectrum of CNCbl 

The photoabsorption spectrum of CNCbl is characterized 
by several distinctive peaks^. The lowest-energy peak in the 
spectrum occurs at « 2.25 eV and it is associated with the so- 
called electronically excited state of CNCbl. The largest 
peak, on the other hand, is located at ^ 3.4 eV. There have 
been various DFT based studies to identify the origin 

of the peaks in the spectra of the Cbl cofactors. 

For a long time, the lowest-energy part of the photoabsorp¬ 
tion spectrum of Cbl had been associated with a transition 
which is dominated by the ;r ^ ;r* excitations of the cor- 
rin ring^ which represent electronic transitions among the 
2p7t orbitals in the corrin ring. The reason for this assign¬ 
ment is that the metal-free corrinoids also exhibit a simi¬ 
lar peak in the spectrum. The time-dependent DFT calcu¬ 
lations with the BP86 functionals on the truncated molecule 
Im-[Co^^^(corrin)]-CN+ have supported the view that the 
state corresponds mainly to ;r ^ ;r* transitions in the corrin 
ring^^. On the other hand, the ultrafast transient absorption 
spectroscopy measurements on CNCbl^^’^^ have suggested a 
picture for the transition in which charge transfer from the 
K orbitals of the corrin ring to the Co atom takes place. In 
spite of these developments, today it is considered that much 
still remains to be understood about the nature of the excited 
states of Cbl cofactors. In particular, there is much interest to 
clarify the origin of the lowest excited state. 

In order to compare the photoabsorption spectrum of 
CNCbl with the HF+QMC data on the electronic structure 
of Im-[Co^^^(corrin)]-CN+, we assume that the rigid band 
approximation can be used for describing the single-particle 
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spectral weight distribution of the Haldane-Anderson model. 
According to the HF-i-QMC results shown in Fig. 7 for 
U = 36 eV, when the chemical potential ji ^ —8.5 eV, the 
truncated system Im-[Co^^^(corrin)]-CN+ has 238 electrons. 
In this case, the highest occupied states consist of the Co t 2 g- 
like states which are the 3dxz, 3d^2_y2 and Sdy^ NAO’s in our 
coordinate system. Smaller amount of weight from the 2p7r 
orbitals of the corrin ring is also included. In the dominant 
lowest-energy photoabsorption process, an electron initially 
occupying a Co r 2 g-like state would be transferred to the low¬ 
est excited state located at ^ —5.5 eV. We have seen in Sec¬ 
tion III.C that this state at ^ — 5.5 eV corresponds to an im¬ 
purity bound state and consists mainly of the 2p7r NAO’s at 
the CN axial ligand. The impurity bound state contains lesser 
amounts of spectral weight from the 2p7r NAO’s of the corrin 
ring as well as from the Co(3dxy) NAO. Hence, according to 
HF-fQMC, the lowest excited state in the absorption spectrum 
corresponds mainly to a transition involving the transfer of an 
electron from the Co r 2 g-like states to the impurity bound state, 
which consists of states at the CN axial ligand, the corrin ring 
and a Co eg NAO. 

If the electron is excited to the second lowest energy of Im- 
[Co^^^(corrin)]-CN+, then the in-gap states located at —4 eV 
would become occupied. This impurity bound state con¬ 
sists mainly of the 2pG and 2s NAO’s of the CN axial lig¬ 
and. Smaller amounts of spectral weight from the 2pG or¬ 
bitals of the corrin ring and of the Co(3dy^2_p) NAO are also 
included. Hence, this electronic transition would be domi¬ 
nated by the transfer of an electron from the Co r 2 g-like states 
to the second impurity bound state. The transitions to the 
« —5.5 eV and « —4 eV impurity bound states would be the 
main features of the photoabsorption spectrum expected for 
Im-[Co^^^(corrin)]-CN+ according to the HF-fQMC data for 
the Haldane-Anderson model. 

In summary, according to the HF-fQMC calculations the 
lowest-energy excitations are dominated by electron trans¬ 
fer from the Co r 2 g-like states to the CN axial ligand. Now, 
the HF-fQMC approach might be too simple to correctly de¬ 
duce the excitation spectrum of vitamin B 12 . Nevertheless, 
the HF-fQMC results presented here suggest that it would be 
useful to look for signatures of a possible 3d CN ligand 
electron transfer process in the excitation spectrum of CNCbl. 
According to the HF-fQMC calculations, such 3d —> CN lig¬ 
and transitions would be accompanied by the vanishing of the 
antiferromagnetic correlations between the Co eg and the CN 
ligand NAO’s along with the collapse of the magnetic mo¬ 
ments at the CN ligand. In addition, when an electron is re¬ 
moved from r 2 g-like orbitals, a magnetic moment would de¬ 
velop due to the lifting of the double occupancy in these or¬ 
bitals. Hence, the 3d -^CN ligand electron transfer would be 
accompanied by the magnetic moment transfer in the reverse 
direction from the CN ligand to the Co t 2 g orbitals. It would 
be useful to probe experimentally these magnetic correlations 
predicted by HF-fQMC. In particular, it would be interesting 
to see how the many-body effects discussed here would influ¬ 
ence the magnetic circular dichroism spectrum of CNCbl^^. 


C. DFT+QMC approach 

An alternative method for studying the electronic structure 
and correlations of organometallic molecules such as Cbl is 
the DFT-fQMC approach. There are various differences be¬ 
tween the HF-i-QMC and DFT-i-QMC approaches. For exam¬ 
ple, in HF-fQMC the correlations of the and p orbitals in 
the molecules are treated at the HF level rather than the more 
accurate DFT. In addition, in HF-fQMC a bare unscreened 
Coulomb interaction is used in the QMC part. The bare 
Coulomb repulsion has a value of order 36 eV when evaluated 
in the basis of the atomic orbitals with the 6-3IG Gaussian 
functions. In the DFT-fQMC, the QMC part of the calcula¬ 
tions are performed by using a renormalized Coulomb inter¬ 
action, which is usually of order 4 eV. This renormalization 
is due to screening by the 5' and p electrons in the system, 
and the renormalized Coulomb interaction can be obtained by 
the constrained local-density approximation^"^’^^’^^. It is im¬ 
portant to note that in both approaches the choices used for 
the Coulomb interaction and the double-counting poc repre¬ 
sent different levels of approximations'^^. Hence, it would be 
useful to compare them. 

A transition-metal atom placed in an organic molecule, as 
in Cbl, is one of the simplest strongly-interacting systems. 
With the availability of spectroscopic data, this is a case where 
the combined electronic-structure and many-body approaches 
can be tested. In this respect, a thorough comparison of the 
HF-fQMC and DFT-fQMC approaches would be useful. Be¬ 
cause of manuscript length considerations, this will be pre¬ 
sented elsewhere^^. 


D. Inter-orbital Coulomb interaction and the Hund’s coupling 


In this paper, we have presented HF-i-QMC results, where 
we used a constant U for the intra-orbital Coulomb interac¬ 
tion in the QMC part. It is well known that the inter-orbital 
Coulomb interaction is important in the cobaltates^^. Here, 
the inter-orbital Coulomb interactions are taken into account 
only at the Hartree-Fock level. A more accurate approach 
would be to include the inter-orbital Coulomb interaction in 
the QMC calculations. This can be done by replacing the last 
term in Eq. (1) by 


1 

2 


v,v' ^ 



Jyy') ^VG^y^ G 


( 22 ) 


where Uyy/ and Jyy/ are the Coulomb matrix elements for the 
direct and the Hund’s exchange couplings. For consistency, it 
is necessary to evaluate the Coulomb matrix elements in the 
NAO basis, since the one-electron parameters of the Anderson 
Hamiltonian are obtained by using the NAO’s. 

The interaction, Eq. (22), includes only density-density 
terms, hence it can be implemented in the Hirsch-Eye QMC 
algorithm. This improvement would still neglect the pair¬ 
hopping and the spin-flip terms of the Hund’s coupling. 
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which have been shown to be important in transition-metal 
oxides^ The double-counting correction /Idc would 

also need to be re-evaluated after including the inter-orbital 
Coulomb interactions'^^"'^^. 

We think that it is necessary to include the inter-orbital 
Coulomb interaction before performing quantitative compar¬ 
isons with the experimental data on the Cbl cofactors. The 
Hund’s coupling may influence the magnetic correlations and 
cause additional splitting of the energy levels. After these im¬ 
provements, the spin state of Co can be determined. In addi¬ 
tion, the uniform magnetic susceptibility for electronic spins 
can be calculated to extract the effective magnetic moment for 
the molecule. 


V. SUMMARY AND CONCLUSIONS 

In summary, we have studied the electronic structure and 
correlations of CNCbl from the perspective of the many-body 
physics. For this purpose, we have used the framework of 
the multi-orbital single-impurity Haldane-Anderson model of 
a transition-metal impurity placed in a semiconductor host. 
First, we have constructed an effective Haldane-Anderson 
model by using the HF approximation. In particular, we de¬ 
termined the one-electron parameters of the Anderson Hamil¬ 
tonian from the Fock matrix written in the basis of the natu¬ 
ral atomic orbitals by using HF. The double-counting of the 
Coulomb interaction within HF-fQMC was taken into account 
by a chemical shift /ioc of the Co(3Jv) levels. We used an 
orbital independent intra-orbital Coulomb interaction U in the 
QMC calculations neglecting the inter-orbital terms. We pre¬ 
sented QMC results for various values U. However, we con¬ 
centrated on the case U = 36 eV, since this yields values for 
the HOMO-LUMO gap and the total Co(3 J) electron number 
which are comparable to the experimental data. 

The QMC data on this effective model showed how the en¬ 
ergy gap found in the HF approximation is reduced by the 
generation of new in-gap states. The correlated nature of the 
induced states is clearly seen. In particular, states arising from 
the double-occupancy of the Co t 2 g orbitals are induced near 
the gap edge above the HOMO level from HF, whereas im¬ 
purity bound states are induced below the LUMO level. The 
impurity bound states arise from the strong hybridization of 
the Co Cg orbitals with the host eigenstates and contain sig- 
niflcant amount of weight from the CN axial ligand and the 
corrin ring. We have also seen that magnetic moments can 
develop at the CN axial ligand and in the corrin ring, which 
are antiferromagnetically correlated with the Co eg moments. 
These disappear rapidly with electron Ailing of the impurity 
bound states. 

We have also presented a comparison of the HF-fQMC data 
with the DFT calculations on Im-[Co^^^ (corrin)]-CN+. While 
these two approaches yield similar values for the HOMO- 
LUMO gap and the total Co(3 J) occupation number, the re¬ 
sults on the distribution of the single-particle spectral weights 
can be very different. In particular, the Hubbard-type splitting 
of the Co(3dy) states, the impurity bound states and the mag¬ 
netic correlations are not contained in the DFT results. We 


have also presented the predictions of HF-fQMC for the pho¬ 
toabsorption spectrum according to which the lowest-energy 
photoabsorption excitations are dominated by electron trans¬ 
fer from the Co t 2 g orbitals to the impurity bound states. 

It will be interesting to perform similar calculations for 
AdoCbl and MeCbl and compare their results with the results 
presented here on CNCbl. In particular, it would be useful to 
see whether it is possible to understand the catalytic function¬ 
ing of these Cbl cofactors by using the HF-fQMC calculations. 
It would also be useful to use the same framework to study 
hemoglobin containing Fe. 

It is worth noting that the Co-C bond in Cbl is the first case 
of a metal-carbon bond to be found in biology. All of the 
known reactions of the Cbl-dependent enzymes involve the 
making and breaking of the Co-C bond^. In this respect, it is 
interesting that the impurity bound state contains significant 
amount of spectral weight from the CN axial ligand attached 
to Co. It would also be interesting to see whether the impurity 
bound state plays a role in the important Co-C bonding. If 
the making and breaking of the Co-C bond somehow involves 
the magnetic correlations around the Co atom, then this could 
be a way such that the occupancy of the impurity bound state 
becomes important. 

The present HF-fQMC approach contains various approx¬ 
imations. For instance, the inter-orbital Coulomb interaction 
is not treated by QMC, instead it is included at the Hartree- 
Fock level. In addition, the bare Coulomb repulsion U and 
the double-counting /Idc are used more as variables than as 
true ab initio parameters. The geometrical parameters for 
the molecular structure are simple estimates since the actual 
atomic coordinates for CNCbl are not known. In spite of 
these, the HF-fQMC yields reasonable values for the HOMO- 
LUMO gap and the total Co(3d) occupation number, and in¬ 
teresting results for the photoabsorption spectrum. In partic¬ 
ular, the HF-i-QMC results emphasize the many-body effects 
such as the Hubbard splitting of the Co(3dy) states and the 
possible existence of impurity bound states and magnetic cor¬ 
relations. 

In conclusion, according to the HF-fQMC results presented 
here, there is a possibility that an impurity bound state exists 
in the electronic spectrum of CNCbl. We think that it is neces¬ 
sary to seek further experimental evidence for this prediction 
of the HF-fQMC calculations. 
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